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£f) • Perturbative QCD uses the Faddeev-Popov gauge-fixing procedure, which leads 

C D ' to ghosts and the local BRST invariance of the gauge-fixed perturbative QCD 

£Nj action. In the asymptotic regime, where perturbative QCD is relevant, Gribov 

T— H , copies can be neglected. In the nonperturbative regime, one must adopt either a 

£Nj 1 nonlocal Gribov-copy free gauge (e.g., Laplacian gauge) or attempt to maintain 

C ~ ) ' local BRST invariance at the expense of admitting Gribov copies. These issues 

| are explored and discussed. In addition, the relationship between recent Dyson- 

Schwinger based calculations of the infrared behavior of QCD Green's functions 
and the lattice calculation of these quantities is examined. 

1. Introduction 

> ' 

. i-h . Perturbative quantum chrome-dynamics (QCD) is formulated using the 

' Faddeev-Popov gauge-fixing procedure, which introduces ghost fields and 

leads to the local BRST invariance of the gauge-fixed perturbative QCD 
action. These perturbative gauge fixing schemes include, e.g., the standard 
choices of covariant, Coulomb and axial gauge fixing. These are entirely 
adequate for the purpose of studying perturbative QCD, however, they fail 
in the nonperturbative regime due to the presence of Gribov copies. Per- 
turbative QCD works because in doing a weak-field expansion around the 

= configuration these Gribov copies are not encountered . 

One could define nonperturbative QCD by imposing a non-local Gribov- 
copy free gauge fixing (such as Laplacian gauge) or, alternatively, one could 
attempt to maintain local BRST invariance at the cost of admitting Gribov 
copies. One of the well-known difficulties for the latter option is the prob- 
lem of pairs of Gribov copies with opposite sign giving a vanishing path 
integral 2 ' 3 ' 4 - 5 . Whether or not a local BRST invariance for QCD can be 
maintained in the nonperturbative regime remains an open problem. 

The standard lattice definition of QCD is equivalent to the choice of a 
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Gribov copy free gauge-fixing. There is a negligible chance of selecting two 
gauge-equivalent configurations (strictly zero except for numerical round- 
off error). Calculations of physical observables are unaffected by arbitrary 
gauge transformations on the configurations in the ideal gauge-fixed en- 
semble. A lattice QCD calculation using an ideal gauge-fixed ensemble will 
give a result for a gauge- invariant (i.e., physical) quantity which is identi- 
cal to doing no gauge fixing at all, i.e., equivalent to the standard lattice 
calculation of physical quantities. 

We begin by reviewing the standard arguments for constructing QCD 
perturbation theory, which use the Faddeev-Popov gauge fixing procedure 
to construct the pcrturbative QCD gauge-fixed Lagrangian density. The 
naive Lagrangian density of QCD is £qcd = —jF^F^ + J2f Qf{ip ~ 
rrif)qf, where the index / corresponds to the quark flavours. The naive 
Lagrangian is neither gauge-fixed nor renormalized, however it is invariant 
under local SU(3) C gauge transformations g(x). For arbitrary, small ui a {x) 
we have g(x) = aq>{-ig s (A a /2) u a {x)} G SU(3), where the A a /2 = t a are 
the generators of the gauge group SU(3) and the index a runs over the 
eight generator labels a = 1, 2, 8. 

Consider some gauge-invariant Green's function (for the time be- 
ing we shall concern ourselves only with gluons) (Q\ T(0[A]) = 
J DA 0[A] e lS ^/ J DA e lS ^ , where 0[A] is some gauge-independent 
quantity depending on the gauge field, A^x). We see that the gauge- 
independence of 0[A] and S[A] gives rise to an infinite quantity in both the 
numerator and denominator, which must be eliminated by gauge-fixing. 
The Minkowski-spacc Green's functions arc defined as the Wick-rotated 
versions of the Euclidean ones. 

The gauge orbit for some configuration A^ is defined to be the set of all 
of its gauge-equivalent configurations. Each point A^ on the gauge orbit is 
obtained by acting upon A^ with the gauge transformation g. By definition 
the action, S[A], is gauge invariant and so all configurations on the gauge 
orbit have the same action, e.g., see the illustration in Fig. 1. 

2. Gribov Copies and the Faddeev-Popov Determinant 

Any gauge-fixing procedure defines a surface in gauge-field configuration 
space. Fig. 2 is a depiction of these surfaces represented as dashed lines 
intersecting the gauge orbits within this configuration space. Of course, in 
general, the gauge orbits are hypersurfaces as are the gauge-fixing surfaces. 
Any gauge-fixing surface must, by definition, only intersect the gauge orbits 
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orbit 



Figure 1. Illustration of the gauge orbit containing and indicating the effect of 
acting on A M with the gauge transformation g. The action S[A] is constant around the 
orbit. 

at distinct isolated points in field configuration space. For this reason, it 
is sufficient to use lines for the simple illustration of the concepts here. An 
ideal (or complete) gauge-fixing condition, F[A] = 0, defines a surface called 
the Fundamental Modular Region (FMR) that intersects each gauge orbit 
once and only once and typically where possible contains the trivial config- 
uration = 0. A non-ideal gauge-fixing condition, F'[A] = 0, defines a 
surface or surfaces which intersect the gauge orbit more than once. These 
multiple intersections of the non-ideal gauge fixing surface(s) with the gauge 
orbit are referred to as Gribov copies 2 - 3 ' 4 - 5 Lorentz gauge (d fl A fl (x) = 0) 
for example, has many Gribov copies per gauge orbit. By definition an 
ideal gauge fixing is free from Gribov copies. The ideal gauge-fixing surface 
F[A] = specifies the FMR for that gauge choice. Typically the gauge 
fixing condition depends on a space-time coordinate, (e.g., Lorentz gauge, 
axial gauge, etc.), and so we write the gauge fixing condition more generally 
as F([A];x) = 0. 




F'[AJ=0 



Figure 2. Ideal, and non-ideal, F'[A], gauge-fixing. 

Let us denote one arbitrary gauge configuration per gauge orbit, A^, as 
the origin for that gauge orbit, i.e., corresponding to g = on that orbit. 
Then each gauge orbit can be labelled by A^ and the set of all such A^ is 
equivalent to one particular, complete specification of the gauge. Under a 
gauge transformation, g, we move from the origin of the gauge orbit to the 
configuration, A 9 ^, where by definition A° — ► A^ = gA^g^ — (i/g s )(d ll g)g ji ■ 
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Let us denote for each gauge orbit the gauge transformation, g = g[A°], 
as the transformation which takes us from the origin of that orbit, A^, 
to the corresponding configuration on the FMR, yl FMR = A 9 , which is 
specified by the ideal gauge fixing condition x) — 0. In other words, 

an ideal gauge fixing has a unique g which satisfies F([A 9 ];x)\g — and 
hence specifies the FMR as A' 9 = A FMR G FMR. Note then that we have 
fVA = JVA°J Vg = J VA™ R J V(g g) . 

The inverse Faddeev- Popov determinant is defined as the integral over 
the gauge group of the gauge-fixing condition, i.e., 



A^L4 FMR ] ee J Vg 6[F[A\] = J Vg 5(g ~g) 



det 



( SF({A];x y 
\ Sg(y) 



(1) 



Let us define the matrix M F [A] as M F ([A];x,y) ab = 6F a ([A];x)/6g b (y) . 
Then the Faddeev-Popov determinant for an arbitrary configuration A^ 
can be defined as A;?L4] = |detMpL4]|. (The reason for the name is now 
clear). Note that we have consistency, since A^ 1 [A FMR ] = A^ 1 ^] = 
JVg 8{g-~g)A F 1 [A\. 

We have 1 = j Vg A F [A] S[F[A}} = JV(g- g) A F [A] S[F[A}} by defi- 
nition and hence 

J VA FMR = J VA FMR J V(g - g) A F [A]8[F[A]]= J VA A F [A]6[F[A\] (2) 

Since for an ideal gauge- fixing there is one and only one g per gauge orbit, 
such that F(L4];x)|g = 0, then detM^^] is non-zero on the FMR. It 
follows that since there is at least one smooth path between any two con- 
figurations in the FMR and since the determinant cannot be zero on the 
FMR, then it cannot change sign on the FMR. The first Gribov horizon is 
defined to be those configurations with det M F [A] — which lie closest to 
the FMR. By definition the determinant can change sign on or outside this 
horizon. Clearly, the FMR is contained within the first Gribov horizon and 
for an ideal gauge fixing, since the sign of the determinant cannot change, 
we can replace | det M F \ with det M F , [i.e., the overall sign of the functional 
integral is normalized away in the ratio of functional integrals]. 

These results are generalizations of results from ordinary calculus, where 
|det {dfi/dxj)\~~^ a = j dxi ■ ■ ■ dx n 5^ n \j(x)) and if there is one and only 

one x which is a solution of ,f{x) = then the matrix Mij = dfi/dxj 
is invertible (i.e., non-singular) on the hypersurface f(x) — and hence 
det M ^ 0. 
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3. Generalized Faddeev-Popov Technique 

Let us now assume that we have a family of ideal gauge fixings F([A};x) = 
f([A];x) — c(x) for any Lorentz scalar c(x) and for /(L4];x) being some 
Lorentz scalar function, (e.g., d^A^(x) or n^A^(x) or similar or any non- 
local generalizations of these). Therefore, using the fact that we remain 
in the FMR and can drop the modulus on the determinant, we have 
JVA FMR = JVA detM F L4] 5[f[A] - c] . Since c{x) is an arbitrary func- 
tion, we can define a new "gauge" as the Gaussian weighted average over 
c(x), i.e., 

J VA FMR cx J Vc ex p{-^ J d 4 xc{x) 2 ^ J VA detMp [A] 6[f[A] - c] 
cx J VA detM F [A]exp|-^y d 4 xf([A];x) 2 ^ 
cx J VAV X V X cxp j-i J d 4 xd 4 y x{x)M F ([A];x,y) X (y)} 

xexp|-^|dV([^]^) 2 } , (3) 

where we have introduced the anti-commuting ghost fields \ and \. Note 
that this kind of ideal gauge fixing does not choose just one gauge config- 
uration on the gauge orbit, but rather is some Gaussian weighted average 
over gauge fields on the gauge orbit. We then obtain 

where 

-\F^F« V - 1 (f([A];x)f + J2q f (ip- m f )q f 



Si[q,q,A,x,x] = ^ d * x 



+ J d^xdS x(x)M F ([A};x,y) X (y) . (5) 



4. Standard Gauge Fixing 

We can now recover standard gauge fixing schemes as special cases of this 
generalized form. First consider standard covariant gauge, which we obtain 
by taking /(L4];a;) = d ll A lx {x) and by neglecting the fact that this leads 
to Gribov copies. We need to evaluate Mf[^4] in the vicinity of the gauge- 
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fixing surface (specified by g): 

Under an infinitesimal gauge transformation about the FMR, Sg = g — g, 
we have -> (A^^)^ where 

(A* +s X(x) = (AX(x) + g s f abc w b {x)A^{x) d^ a {x) + 0(c 2 ) (7) 

and hence in the neighbourhood of the gauge fixing surface (i.e., for small 
fluctations along the gauge orbit around A™ R ), we have 



M F ([A];x,yy b =d* 6Aa " {x) 



(8) 

w=0 



^6(Su b (y)] 

= d£ ( [-d^5 ab + g s r bc A^(x)}8^(x - y)) . 
We then recover the standard covariant gauge-fixed form of the QCD action 



Sz[q,q,A,x,x] = J d * x 



+(d»Xa){d»5 ab -gfabcA^Xb- (9) 

However, this gauge fixing has not removed the Gribov copies and so the 
formal manipulations which lead to this action are not valid. This Lorentz 
covariant set of naive gauges corresponds to a Gaussian weighted average 
over generalized Lorentz gauges, where the gauge parameter £ is the width 
of the Gaussian distribution over the configurations on the gauge orbit. 
Setting £ = we see that the width vanishes and we obtain Landau gauge 
(equivalent to Lorentz gauge, d^A^(x) — 0). Choosing £ = 1 is referred to 
as "Feynman gauge" and so on. We can similarly derive the QCD action 
for axial gauge. 



5. Discussion and Conclusions 

There is no known Gribov-copy-free gauge fixing which is a local function 
of A^(x). In other words, such a gauge fixing cannot be expressed as a 
function of A^(x) and a finite number of its derivatives, i.e., F(L4];a;) ^ 
F(d )Jl , A fl (x)) for all x. Hence, the ideal gauge-fixed action, in 
Eq. (5) becomes non-local and gives rise to a nonlocal quantum field theory. 
Since this action serves as the basis for the proof of the renormalizability 
of QCD, the proof of asymptotic freedom, local BRST invariance, and the 
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Schwinger-Dyson equations 6,7 (to name but a few) the nonlocality of the 
action leaves us without a first-principles proof of these features of QCD in 
the nonperturbative context. 

The lattice implementation of Landau gauge finds local minima of the 
gauge fixing functional, which correspond to configurations lying inside the 
first Gribov horizon. The remaining Gribov copies after this partial gauge 
fixing then necessarily all have the same sign (positive) for the Faddeev- 
Popov determinant and hence add coherently in the functional integral. 
This ensures that the ghost propagator is positive definite 7,8 . The deriva- 
tion of the Dyson-Schwinger equations is based on the fact that the in- 
tegral of a total derivative vanishes 6 provided that the surface integral of 
the integrand vanishes when integrated over the boundary of the region. 
Since the Faddeev-Popov determinant vanishes on the first Gribov horizon, 
then we can still derive the standard Dyson-Schwinger equations from the 
Landau gauge fixed QCD action even if we restrict the gauge fields to lie 
within the first Gribov horizon. This is equivalent to requiring that the 
ghost propagator be positive definite. Thus it is valid to compare lattice 
Landau-gauge calculations with Dyson-Schwinger based calculations (with 
a positive definite ghost propagator), since these both consist of consider- 
ing configurations within the first Gribov horizon. An extensive body of 
lattice calculations exist for the Landau gauge gluon 9 and quark 10,12 prop- 
agators and most recently for the quark- gluon vertex 13 . Similarly, calcula- 
tions in Laplacian gauge (an ideal gauge) fixing have also recently become 
available 14,15 . 

It is well-established that QCD is asymptotically free, i.e., it is weak- 
coupling at large momenta. In the weak coupling limit the functional 
integral is dominated by small action configurations. As a consequence, 
momentum-space Green's functions at large momenta will receive their 
dominant contributions in the path integral from configurations near the 
trivial gauge orbit, i.e., the orbit containing = 0, since this orbit min- 
imizes the action. If we use standard lattice gauge fixing, which neglects 
the fact that Gribov copies are present, then at large momenta J DA will 
be dominated by configurations lying on the gauge-fixed surfaces in the 
neighbourhood of each of the Gribov copies on the trivial orbit. Since for 
small field fluctuations the Gribov copies cannot be aware of each other, we 
merely overcount the contribution by a factor equal to the number of copies 
on the trivial orbit. This overcounting is normalized away in the ratio of 
functional integrals. Thus it is possible to understand why Gribov copies 
can be neglected at large momenta and why it is sufficient to use standard 
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gauge fixing schemes as the basis for calculations in perturbative QCD. 
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